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Abstract 

We extract transition amplitudes among matter constituents of the universe from the 
solutions of the Wheeler De Witt equation. The physical interpretation of these solutions 
is then reached by an analysis of the properties of the transition amplitudes. The 
interpretation so obtained is based on the current carried by these solutions and confirms 
ideas put forward by Vilenkin. 
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1 Introduction 



In the hamiltonian formulation of General Relativity coupled to matter fields, the Euler- 
Lagrange equations split into 6 dynamical equations plus the equations for the matter 
fields and 4 constraints which express the invariance of the theory under local changes of 
coordinates. Therefore, when quantizing this system in a gauge invariant manner, i.e. a 
la Dirac, the "wave function" must be annihilated by these four constraints that are now 
operator valued. This clear structure dictated by the reparametrization invariance turns 
out to lead to very complicated problems concerning the interpretation of this "wave 
function". By interpretation, one designates the problem of extracting (probabilistic) 
predictions concerning the evolution of matter and gravity configurations. The origin 
of the difficulties can be blamed on the absence of time and the accompanying unitary 
evolution on which quantum mechanics is based[l]H[PJP|H]. 

In order to describe more precisely some of these difficulties, from now on we shall 
pursue the discussion in mini-superspace. This drastic restriction to homogeneous and 
isotropic three-geometries offers the double advantage of removing the U.V. problem 
that plagues the local theory while keeping the problem concerning the interpretation 
of the wave function. In this restricted configuration space, gravity is described by the 
scale factor a, matter by homogeneous fields that we shall denoted collectively by and 
the wave \l/(a, 0) is constrained to satisfy a single global Wheeler De Witt (WDW) equa- 
tion. The questions are the following: How to proceed to read from \l/(a, 0) predictions 
concerning quantum events such as, for instance, transitions rates among the fields ? 
and: Are those transitions described by a unitary evolution ? 

Both questions have received attention and many different schemes have been proposed M. 
Let us mention here only the scheme based on the hypothesis that \^f(a, <f))\ 2 does posses 
a probabilistic interpretation || (at least in terms of conditional probabilities!]) and the 

scheme based on the conserved current ty*id$. The answers to the second question 
are partially related with this choice and range from, "Yes, the evolution is unitary 0" , 
to "no, unitarity is violated [f|]", and includes the middle attitude: "unitarity is only ap- 
proximatively conserved ||". The peculiar aspect of these widely distributed answers is 
that they arise from the same starting point: the WDW equation and its solutions. The 
disagreements build up with the choice of the treatment required to extract information 
from \l/(a, 0). 

In the present article, we clarify the mathematical aspects of these treatments by 
analysing matter interactions from the solutions of the WDW equation. More precisely, 
by studying transitions, we identify the coefficient C n (a) that weights the n-th state 
at a and that replaces the amplitude c n (t) to be in the n-th matter state at time t 
in conventional quantum mechanics. The unambiguous identification is based on two 
criteria: 1) In the absence of transitions, C n (a) must be constant. 2) When one simplifies 
the equation governing their dependence in a by treating gravity in the background field 
approximation (BFA), the resulting equation must be the Schroedinger equation. This 
mathematical procedure of extracting the C n (a) is based on || || [|l(| and does not require 
an a priori interpretation of \I/(a, 0). On the contrary, our program is first to determine 
the properties of the coefficients C n (a) and only then to examine the question of its 
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interpretation in the light of these properties. 

They reveal the existence of three regimes which are delineated by the values of the 
parameters describing matter transitions in quantum cosmology. In the case of weak 
interactions occurring close to equilibrium in a macroscopic ||fll], the coefficients C n (a) 



are equal to the amplitudes c n (t), solutions of the corresponding Schroedinger equation 
evaluated in the geometry described by a(t). In the second regime, the departure from 
equilibrium and/or the importance of the interactions lead to C n (a) that no longer co- 
incides with c n (t). However when gravity is still correctly described by WKB waves, 
the C n (a) still satisfy the "unitary" equation J2 n l^n(a)| 2 = 1, up to negligeable correc- 
tions. In the third regime, the interactions are so violent that the propagation of a is 
affected by the quantum transition acts. In that case, J2 n \Cn(a)\ 2 1- This "violation" 
is a direct manifestation of the modification of the propagation of a by the transitions 
themselves. It is also kinematically related to the conservation of the current carried 
by \l/(a,0). Under these extreme conditions, there is no possibility of defining a back- 
ground. Neither, therefore, should there be any possibility of interpreting \l/(a,0) using 
the conventional rules of quantum mechanics. This does not mean that no predictions 
can be made, it simply means that the conventional analysis cannot be performed in 
situ. By propagating $?(a, (ft) outside this regime, one can then perfectly determine its 
physical outcome. 

In conclusion, in this article, we shall show that when one requires that the con- 
ventional description of matter transitions is recovered from quantum cosmology, the 
probabilistic interpretation of ^(a, <p) must be based on its current and not on its norm. 
We shall also show that the interpretation of C n (a) as the amplitude of probability to 
find the n-th matter state at a is valid as long as the propagation of gravity is not sig- 
nificantly affected by the interactions. Both aspects have been put forward by Vilenkin 
in [||]. However, to the knowledge of the author, they have never been made as explicit 
as in the present paper. Furthermore, in contradistinction to ||, our small parameter is 
the coupling constant among the quantum systems and not their energy. This allows to 
reach more general conclusions. 



2 The identification, the evolution and the meaning 
of the coefficients C n (a) 

As said in the Introduction, we shall use perturbation theory applied to matter inter- 
actions as a guide to identify the coefficients C n (a). Before accomplishing this program, 
we briefly present the kinematical properties at work in Quantum Cosmology, in the 
absence of these interactions, see || for more details. 

The coefficients in absence of interactions 

For simplicity, the matter system is chosen in such a way that the free hamiltonian H^a) 
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does not lead to transitions. Thus, the general solution of the Schroedinger equation is 

xM) = £c n e-^'^(O (0 | n) (1) 

n 

where the eigenstates of the free hamiltonian H^(a) satisfy 

H° m (a(t))\n) = E n (a(t))\n) (2) 

Their time dependence arises only through the equation governing the background prop- 
agation a = a(t). 

The coefficients c n are interpreted as the amplitudes of probability to find the mat- 
ter system in the n-th state. Therefore, by convention, they are normalized so that 
En l c n| 2 = 1- Furthermore, in the present case, they are constant. In this respect, notice 
that one must consider either interactions with the external world, or self interactions, 
or non adiabatic transitions [0] in order to give physical substance to the probabilistic 



interpretation of c n . Indeed, in the absence of interactions, no interference among the 
c n will show up. 

We shall now determine to what extend these properties are recovered in quantum 
cosmology from the solutions of the Wheeler De Witt equation. 

In minisuperspace, when matter is characterized by an energy E n (a), the gravita- 
tional part of the action satisfies the Hamilton- Jacobi constraint 

H G (a) + E n (a) = h E n (a) = (3) 

where G is Newton's constant, k is equal to ±1 or for respectively open, closed and 
flat three surfaces and A is the cosmo logical constant. The solution of this equation is 
simply S n (a) = J a da'p n (a') where the momentum of a driven by the E n (a) 



p n (a) = tG-^ko? + Aa 4 + 2GaE n (a) (4) 

where =F correspond respectively to expanding and contracting universes. 

Upon quantizing gravity, the Hamilton- Jacobi constraint becomes the operator val- 
ued equation (the WDW equation) 

# G (a)+i£(a>)]*(a,0) = O (5) 

When the quantum matter states are characterized by constants of motion, \l/(a, 0) can 
be decomposed as 

^(a, ( p) = Y / Cn^(a;n)( ( p\n) (6) 

n 

where the n-th gravitational wave, entangled^] to the n-th matter state, is a solution of 

G 2 d 2 a + kcl 2 + Aa 4 + 2GaE n (a)} (a; n) = (7) 



1 This is the main reason for which we have chosen to work with matter states characterized by 
constants of motion. It allows an unambiguous decomposition of the total wave in terms of products. 
This is to be compared with the difficulties to perform this decomposition in the general case wherein 
no clear principle seems to exist, see e.g. J| after eq. (2.35) "We choose D in such a way that the 
equations become simple" . Very important also is the fact that our decomposition keeps the linearity 
of the WDW equation when used in a perturbative treatment. 
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This equation is second order in d a and has therefore two independent solutions. This 
has to be the case since classically we can work either with expanding or contracting 
universes. Indeed when using the WKB approximation 



e * J a p n (a')da' 

*(a;n)= (8) 
p\Pn(a)\ 

one verifies that solutions with positive (negative) wronskian 

W n = **{a-n) S tt(a;n) (9) 

correspond to expanding (contracting) universes. 

At this stage, several remarks should be made. Firstly, the decomposition of the 
general solution of the WDW equation is performed by using the same set of quantum 
numbers n that the one used in eq. @ wherein gravity was treated in the background 
field approximation (BFA). The enlargement of the dynamics to a is compensated by 
the WDW constraint. Secondly, it is now through the sign of the Wronskian rather than 
at the classical level that one now chooses expanding or contracting universes. Thirdly, 
with our definition of \I/(a;n), the coefficients C n are constant. 

Having recall these kinematical properties, we can now formulate precisely our math- 
ematical claim |J: in order for C n to satisfy the Schroedinger equation when gravity is 
treated in the BFA, all Wronskians W n must be equal to the same constant. This can 
be already guessed by considering the conserved current carried by \&(a, 4>) 

J= fdcf> V*(a, 0)zU*(a, 0) = E \ C n\ 2 W n (10) 

To work with both J = 1 and unit Wronskians suggests an identification of C n with c n . 
However, at this moment, no physically relevant conclusion?] can be made concerning a 
probabilistic interpretation of the C n . Indeed it is mandatory to consider self interactions 
leading to transitions in order for the C n to vary and interfere, see the remark made after 
eq. (@). 



The a-dependence of the C n 

We return for a moment to quantum mechanics and consider the time dependent per- 
turbation theory for allowing a comparison with the corresponding equations derived in 
quantum cosmology. Upon introducing an interacting hamiltonian H int possessing non 

2 It is nevertheless interesting to compare the present treatment based on the current J to the one in 
which it is the norm that determines probabilities. In that latter case, by definition, the probability to 
be in the n-th state is given by | (n|*(a, <j>)} \ 2 /(^\^} = |C n | 2 p~ 1 (a)/J^ m \C m \ 2 /p^(a), where we have 
used the WKB form for the waves ^(a;n). This "probability" depends on a through the a-dependent 
norm of each ^(a;n), a feature that I find unattractive. Notice however, that in the doubled limit of 
well grouped states living in a macroscopic universe, this dependence in a vanishes, see latter in the text 
for more precision concerning these limits. 
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vanishing matrix elements (n\H int \m) , the time dependence of the coefficients c n is given 
by 

id tCn {t) = £ (n\H mt \m) c m (t) e -*r*!\B~V)-B»V)\ 

m 

Since this equation is linear in c n and first order in id t , when the hamiltonian H int is 
hermitian, one obtains S n |c n (t)| 2 = 1, a necessary condition to keep the probabilistic 
interpretation of |c n (t)| 2 . 

In quantum cosmology, the fact of taking into account the hamiltonian Hi nt is quite 
different from what we just did in usual quantum mechanics wherein one has an external 
time parameter at our disposal. Indeed, Hi nt modifies the propagation of both gravity 
and matter through the modified WDW equation given by 

H G (a) + H° m {a, 0) + H int (a, 0)] E(a, 0) = (12) 

By decomposing the interacting wave S(a, 0) in terms of the free components \l/(a; n) 

S(a,0) = £C m (a) *(a;m)(0|m) (13) 

m 

and by projecting eq. ([TJ) into the bra (n\, we obtain^ 

da$f(a; n) n . , . 1 „ , . a ^ „ , . , . TT . . ^(a: m) . . 

i^a„c„( o) + -gtCM + -£C m (a) <»l^.l-)^f = o (14) 

Since this equation is second order in d a , some analysis is required in order to reveal the 
properties of the evolution it encodes. To this end, we shall first simplify it by making 
use of three approximations. In the next subsection, we shall evaluate the errors they 
induce on the basis of the analysis of . 

The first approximation consists in using the WKB approximation for \l/(a;n). Its 
validity requires d a p n (a) / p^(a) <iC 1, an inequality which is satisfied when the second 



condition is met. This second condition requires that the universe be macroscopic p||TT 
i.e. that the matter sources driving gravity must be macroscopic. By denoting M the 
rest mass of the atoms and Am the energy change induced by the transitions engendered 
by Hi n t, the macroscopic limit guarantees that Am/E n ~ Aiti/NmM <k 1 since Nm, 
the total number of atoms, satisfies Nm 1. Thirdly, we require that the dimensionless 
constant g 2 that characterizes the transition rates be smaller or comparable to unity (g 
is related to H int by (n\H int \m) ~ gAm). 

By applying these three approximations to eq. (JH|), one can drop the second term 
and write the two other terms as 

l d a C n (a) = £ 1 H^|m) Cm{a) elfda/[pm{a ^ Ma/)] M; (15) 



3 Notice that this development does not coincide with the Born-Oppenheimer treatment. What would 
be closer to that treatment, would consist in working with states which diagonalise the total hamiltonian 
H^ n (a,(f) + Hi n t(a,(f)) at fixed a. Together with S. Massar, we shall present this adiabatic treatment 
applied to quantum cosmology in p2] . 
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This equation deserves a few comments. 



Firstly, as eqs. (O, 14), it is linear in C n (a). We point out this fact since many 



approximation schemes P| P |7j] [PH , as the semi-classical treatment, destroy the linearity 
of the WDW equation and therefore the superposition principle. The loss of linearity in 
these treatments results from a quantum averaging performed at an earlier stage. 

Secondly, in order for eq. (p75|) to coincide with eq. (]XTj) in the background field 
approximation, the Wronskians must satisfy W m /W n = 1 for all m,n. The validity of 
the further simplification which consists in treating eq. ( fL5|) in the BFA, requires that 
the C n (a) 7^ be grouped together such that their spread around the mean energy E n 



satisfies (E n — E n )/E n <<c l[Hj|LT]]. Only then can one correctly describe the evolution 



in terms of a single time parameter fl3|1 defined, from the propagation of a, by 



tM = I da 'G^M (16) 



and develop the a-dependent phase of eq. (jig) to first order in E m — E n around E n , see 
|§. Using eq. (|J) and t n (a), one finds, for an expanding universe, 

- / da'(p m (a') - p n (a')) = / dt'(E m (t') - E n (t')) + 0((E m - E n )/E n ) (17) 

Thirdly, when W m /W n = 1, eq. ( |H)| ) leads to "unitary" evolution in the sense 
that E n |C n (a)| 2 = 1. We emphasize that this does not imply that, starting with C n that 
coincide with c n at ao (t = 0), C n (a) will evolve like c n (t(a)). Indeed, as shown in eq. (|17|) , 
the phases of eq. (|15D equal the corresponding phases of eq. ( |TTD only when developed 
to first order in AE and evaluated for the mean energy E n . Therefore, the non-linear 
terms will induce increasing additional phase shifts. (This is not particular to quantum 
cosmology. Indeed, whenever the quantum dynamics is enlarged to a variable formerly 
treated classically, non-linear phase shifts appear, see Then, after a certain lapse 



of t(a), the interferences amongst the C n (a) will posses no relation to those amongst the 
c n (t(a)). Moreover, remote configurations evolve with their own timej5]], see eq. ([11]) for 
the dependence on E n in t n . 

Fourthly, eq. flTo] ) might have been obtained by "first quantizing" the wave S(a, 
In that framework, one postulates that the fundamental equation is 



- z<9 a 2(a,0) = yjH Q {a,4>) + H int (a,(j)) E(a,<p) (18) 

instead of eq. (0), see [I]. Then, by developing the square root to first order in H int 
one obtains eq. ( |15D exactly like eq. (11) is obtained|J. The main weakness of this ad 
hoc approach is that there is no a priori justification to eliminate half of the solutions 
of the Hamilton Jacobi equation before quantization. Furthermore, whether or not it 
offers a good approximation, the importance of the neglected terms cannot be estimated 
without considering the solutions to eq. ([12]). 

Finally, the interpretation of the wave S(a, 0) as the conditional amplitude of prob- 
ability to find ^ at a is rejected by the present analysis of the solutions of eq. (p"2|). 
Indeed, by adopting this interpretation, one would obtain a non-linear equation for the 
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conditional amplitudes since these are non linearly related to C n (a) -recall the presence 
of the normalisation factor p" 1 / 2 (a) stemming from current conservation, see the second 
footnote-. Notice that this would not have been the case if one would have worked with 
the solutions of eq. fll8|) . Notice also that, in contrast to what is presented in ||, our 
derivation of eqs. (|nj, |HS| ) which encode the correlations between matter and gravity 
in no way requires that E(a, 0) be peaked around certain configurations. But it does 
require that gravity be modified by the transitions, otherwise a could not be used to 
parametrize their amplitudes. 

We now address the problem of the corrections to eq. (|l5l). 

The a-dependence of S n |C n (a)| 2 

Our aim is to determine how the terms neglected in passing from eq. ([12]) to eq. (|T5|) 
affect E n |C n (a)| 2 . This is most easily achieved by using the exact conservation law of the 
current carried by S(a,0). Indeed the current carried by E(a, 0) is 



J = J d(j> \E*(a, 4>)iB a E(a, 4>) 

= Y,\ C n\ 2 + J2( C :( a W n (a))\^(a;n)\ 2 = l (19) 



We first notice that there is no systematic departure from unity. Indeed, when the 
interactions cease to act, due for instance to some cooling in an expanding universe, the 
final value of S n |C n (a)| 2 coincides with the initial one. 

Before discussing the reasons that lead to the local departure from unity, we estimate 
the importance of this departure. To this end, it is sufficient to use eq. (|15|) and the WKB 
expression for the free waves \I/(a, <ft). One should also further specify the characteristics 
of the solution under examination. As an example, in the case of a matter dominated 
universe composed of Nm "heavy" atoms of mass M having transitions between inner 
states characterized by energy gaps Am, this correction is 

1 /2 

x , / '.<-* „ \ , . Nl9 n-^n-^^ a(n\Hi nt \m) aN Af Am Am 1 . , 



far from a turning point |16[ and at equilibrium. (The incoherence of the transitions close 
or at equilibrium brings the factor N\[ 2 in the second approximation.) The departure 
from unity is therefore small for three independent reasons: the weakness of the transi- 
tions rates controlled by g, the smallness of the Am/M which allows to neglect the recoil 
of the atoms induced by the transitions and the macroscopic character of the universe 
Nm » I- 

The mechanism that leads to the departure from unity is clear: when the C n de- 
pend on a, they carry some current and therefore their norm should vary accordingly. 
Indeed, when the "potential term" of a second order differential equation varies, the 
induced modifications of the norm and the phase of the solutions are correlated since 
the Wronskian is conserved. The physical origin of these correlated modifications comes 
from the fact that the WDW imposes that the kinetic energy of gravity be modified 



8 



by the presence of H int , see eq. (]T2|). This in turn modifies the norm of the solution. 
Therefore, S ri |C n (a)| 2 7^ 1 is a manifestation of the modification of the propagation of 
gravity induced by the interactions themselves. 

To obtain explicitly the relation between the changes in phase and in amplitude, 
it is instructive to consider the simple case wherein the interactions induce a diagonal 
level shift E n (a) — > E n (a) + 5 n (a). Using eq. (|) in eq. (§), one immediately gets the 
desired relation through the change of the p n (a). Upon considering matter transitions, 
this modification of p n (a) has the following physical meaning. The transition probability 
is given by an integral over a whose norm is ada/p n (a) = dt n , see |§. Therefore, the 
modification of p n (a) is necessary in that it guarantees that the Golden Rule is obtained, 
i.e. that the transition probability increases linearly with the proper time lapse, eq. (|l6l) . 
evaluated in the universe wherein the transition occurs. 

In view of this, it is inviting to work with eigenstates of the total matter hamiltonian 
Hq + Hi n t rather than to work in perturbation theory with the free states. Indeed, 
when working with these eigenstates, one automatically includes the backreaction of the 
(adiabatic part of the) interaction hamiltonian in the definition of the WKB momentum 
p n (a). Since the induced modification of p n (a) will no longer be found, its contribution to 
£ n |C n (a)| 2 7^ 1 will be also absent. This strongly suggests not to interpret £ n |C n (a)| 2 7^ 1 
as a violation of unitarity since the magnitude of the departure from unity depends on 
the perturbative scheme adopted. Together with S. Massar, we shall return to these 
aspects in [II 



Conclusions and additional remarks 

In resume, our analysis of C n (a) shows that there are three regimes delineated by the val- 
ues of the coupling constant g, the relative transition energy Am/M and the macroscopic 
character of the universe, here controlled by Nm- 

In order to be in the first regime, the initial values C n (a ) 7^ must be grouped 
together so that the "mean" time parameter, eq. flTB|), correctly parametrizes the evolu- 
tion of the whole system. One must also requires that the universe be macroscopic and 
that the transitions be not too violent. In that regime, C n (a) = c n (t n (a)). Therefore, 
C n (a) is the amplitude of probability to find the matter system in the n-th state at a. 

In the second regime, the spread of C n (a ) 7^ and/or the violence of the interactions 
and/or the "smallness" of the universe leads to an evolution that cannot be obtained from 
a Schroedinger equation based on a single time parameter. Nevertheless, when the WKB 
approximation for describing the free evolution of gravity is still valid, one still obtains 
a linear equation for the propagation of the C n (a) that guarantees that S n |C n (a)| 2 = 1. 
Moreover, upon considering a set of neighbouring C n (a) for a sufficiently small amount 
of time, the evolution of this restricted set still obeys a Schroedinger equation. Finally, 
as usually in quantum mechanics, remote configurations do not interfere. Therefore, it is 
still perfectly correct to interpret C n (a) as the amplitude of probability to find the n-th 
state at a. Notice however that the decoherence of the background time, that is the fact 
that remote configurations evolve with their own time, would leave no physical meaning 
to mean values of matter operators in which non-interfering configurations would con- 
tribute. Only summations over restricted sets of interfering configurations make physical 
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sense. In this we differ from 00 since, in these works, physical predictions are based 
on the peaks of S(a, </>) which include summation over all states. 

In the third regime, the importance of the interactions leads to backreaction effects on 
the propagation of gravity such that E n |C n (a)| 2 depends on a. A part of this dependence 
comes from the (adiabatic[]l2"||) dressing energy brought in by the interactions. It reflects 
the fact that the relationship between proper time lapse and the momentum p(a) has 
been modified by this dressing energy. Therefore the deviation from unity it engenders 
should certainly not be interpreted as a violation of unitarity. This deviation can be also 
understood from the necessity of applying a "reduction formula" to the universe's wave 
function in order to obtain transition amplitudes amongst its constituents, see ||. 

Far more difficult to interpret are the consequences of the corrections to the WKB ap- 
proximation. Indeed upon dealing with exact solutions for the propagation of a, initially 
expanding (forward) solutions contain, at other radii, backward waves corresponding to 
contracting universes. Moreover, the conservation of the Wronskian tells us that the cur- 
rent of the forward part has increased. The only way to confront this "Klein paradox" 
seems to proceed to the so called third quantization fj]. 

However, we want to emphasize that by adopting this framework, one has not solve 
the question of extracting transition amplitudes occurring within expanding universes. 
Indeed, additional choices must be specified in order to know how to extract from su- 
perpositions of expanding and contracting universes predictions concerning transitions 
in the expanding sector. Without having further specified these choices and having con- 
sidered combinations of forward and backward solutions, it is overhasty to conclude that 
unitarity will (not) be violated in Quantum Cosmology. 
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